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The production rate and kinematics of photons produced in association with Z bosons are studied 
using 2 fb~^ of pp coUision data coUected at the CoUider Detector at Fermilab. The cross section 
for pp + X (where the leptons I are either muons or electrons with dilepton mass Mu > 

40 GeV/c^ , and where the photon has transverse energy i?J > 7 GeV and is well separated from 
the leptons) is 4.6 ± 0.2 (stat) ± 0.3 (syst) ± 0.3 (lum) pb, which is consistent with standard model 
expectations. We use the photon Et distribution from Z7 events where the Z has decayed to fi'^ ii~ , 
e^e~ , or ui? to set limits on anomalous (non-standard-model) trilinear couplings between photons 
and Z bosons. 

PACS numbers: PACS numbers 14.70.Bh, 14.70.Hp, 13.85.Qk 



The study of production [1] is an impor- 
tant test of the standard model (SM) description 
of gauge boson interactions and provides sensitiv- 
ity to physics beyond the SM [2]. The cross 
section is directly sensitive to the trilinear gauge 
couplings at the ZZ7 and vertices, which 
are predicted to vanish in the SM at tree level [ - 
">]. Physics beyond the SM (e.g. compositeness 
or supersymmetry) could alter the cross section 
and the production kinematics. Z7 production is 
also an important background to searches for new 
physics (e.g. in Gauge-Mediated Supersymmetry 
Breaking models ['■]) and Higgs searches. In this 
paper the production properties of Zj events are 
compared to SM predictions, and limits are set 
on anomalous trilinear couplings. 

We present measurements of Zj production 
from pp collisions at y/s — 1.96 TeV at the Teva- 
tron Collider using data obtained with the Col- 
lider Detector at Fermilab (CDF II). We use two 
methods to identify Zj events. For events where 
the Z decays to charged leptons, we first iden- 
tify the Z decay products and then look for the 
associated photon [7]. We use these events to 
measure the Z7 cross section since we are able 
to identify them with low backgrounds even at 
low photon transverse energy (E^) [s]. Anoma- 
lous (non-SM) Z7 couplings (described below) 
would produce an excess of events with high E^, 
so to set limits on these couplings we also in- 
clude events where the Z boson decays into neu- 
trinos. We identify these by looking for events 



with high E'^ and large missing transverse energy 
(^t) ['^]- An integrated luminosity of 2.0 fb~^ is 
used for Z bosons selected in the Z — > 11^ n~ and 
Z ^ vv decay modes and 1.1 fb^^ is used for the 
Z e+e~ decay mode. 

The CDF detector is described in detail else- 
where [10, 11]. The transverse momenta (px) of 
charged particles are measured by an eight-layer 
siHcon strip detector [12] and a 96-layer drift 
chamber (COT) [ ] inside a 1.4 Tesla magnetic 
field. The COT provides tracking coverage with 
high efficiency for \r]\ < 1. Electromagnetic and 
hadronic calorimeters surround the tracking sys- 
tem. They are segmented in a projective tower 
geometry and measure the energies of charged 
and neutral particles in the central {\r]\ < 1.1) 
and forward (1.1 < \r]\ < 3.6) regions. Each 
calorimeter has an electromagnetic shower profile 
detector positioned at the shower maximum [14]. 
The calorimeters are surrounded by muon drift 
chambers covering jryj < 1. Gas Cherenkov coun- 
ters [] ■"] measure the average number of pp in- 
elastic collisions per beam crossing and thereby 
determine the beam luminosity. The readout de- 
cision is made with a fast three-level trigger sys- 
tem that has high efficiency for selecting the Z 
bosons to be used in the offline analysis. 

For the Zj cross section measurement, we se- 
lect Z7 candidate events by identifying Z — > 
/i^/i^ and Z — > e+e~ and then looking for pho- 
ton candidates with E'^ > 7 GeV. Event selec- 
tion starts with inclusive muon (electron) trig- 
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gers requiring muon px > 18 GeV/c (electron 
Et > 18 GeV). Further selection requires the 
trigger muon (electron) to be isolated and to have 
PT > 20 GeV/c {Et > 20 GeV) and to pass 
standard muon (electron) identification criteria. 
Once we have identified a trigger lepton, we re- 
quire a second oppositely charged lepton of the 
same flavor. For the second muon, we accept 
an isolated track with px > 20 GeV/c, |ry| < 1, 
with a minimum-ionizing energy deposit in the 
calorimeter, regardless of whether or not it is as- 
sociated with a reconstructed track in the muon 
chambers. For the second electron, we accept iso- 
lated electromagnetic clusters with Et > 20 GeV 
and 1 7/ 1 < 2.8. One of the electrons is allowed 
to pass a looser isolation requirement, with up 
to three times the amount of additional energy 
near the electron compared to the standard iso- 
lation requirement. The shape of the shower in 
the shower max detector is required to match the 
expected shape for an electron, and in the central 
region we require a match between this shower 
and a track. Since the tracking efficiency falls off 
in the forward region {\ri\ > 1.2), we do not re- 
quire forward electrons to be matched to a track 
and therefore we do not require opposite charge. 
The dilepton invariant mass for both /i'^/i^ and 
e+e^ is required to be larger than 40 GeV/c^. 

We measure trigger and reconstruction efficien- 
cies by studying leptonic decays of W and Z 
bosons. We measure muon trigger efficiencies 
using Z — > i-i^ decays where the two muons 
arc fiducial to different muon detector systems. 
The muon triggers varied during the running pe- 
riod; most of the data were taken with triggers 
that only required a track match in the r-(j) plane 
while a 3-dimensional track match was required 
later during high-luminosity running. The aver- 
age trigger efficiency for muons is (91.6 ± 1.0)%. 
We also use the two muons from Z — >■ /J."*"^" de- 
cays to measure the average muon reconstruction 
and identification efficiency to be (81.3 ± 1.1)%. 

To find the trigger efficiency for a central elec- 
tron that satisfies the selection criteria, we use 
W boson candidates triggered by an electromag- 
netic cluster and but with no tracking re- 
quirement. The trigger efficiency for electrons 
passing the final selection is nearly 100% for elec- 
trons with a high-pT track (> 80 GeV/c), aver- 
ages (95.9 ± 0.6)% for electrons with track pT 
between 20 and 80 GeV/c, and falls sharply for 
electrons with track pT < 14 GeV/ c. We also use 
the two electrons of Z boson candidates to mea- 
sure the electron reconstruction and identifica- 



tion efficiency. This falls slightly during the run 
as the luminosity increases; for central electrons 
it is in the range of 95% to 96% and for forward 
electrons it is in the range of 90% to 92%. 

A small fraction of the Z boson candidates are 
background, primarily from QCD jets that are 
reconstructed as electrons. In the central region 
we use same-sign events to estimate the number 
of background events, which is found to be negli- 
gible in both the muon and the electron datasets. 
In events with one central electron and one for- 
ward electron, the charge of the forward electron 
is not measured reliably, so we normalize to the 
high-mass sideband and find (7.7 ± 1.0)% back- 
ground. 

Once we have selected events with Z bo- 
son candidates, we look for well-identified iso- 
lated [l(i] photons in the central region {\ri\ < 1.1) 
with > 7 GeV that are well separated from 
the Z decay leptons {ARg^ > 0.7). To deter- 
mine the efficiency for identifying real photons, 
we select a sample of Z events where one of the 
final-state leptons has radiated a photon. We 
use tight lepton selection criteria to reduce QCD 
backgrounds and very loose criteria for selecting 
the photon to insure nearly full acceptance for 
photons that have not converted inside the track- 
ing chamber. We reduce the background fraction 
to a negligible level by requiring the three-body 
invariant mass of the two leptons and the pho- 
ton to be 91 ± 5 GeV/c^. Using this method we 
determine that the efficiency for identifying real 
photons is (86 ±2)%, with the largest inefficiency 
coming from the isolation requirement, which is 
needed to reduce QCD backgrounds. 

The largest background to the photons comes 
from QCD jets. A jet can fragment into a high- 
Pt tt^ which decays into two photons. If the two 
photons are sufficiently collinear, their electro- 
magnetic showers in the calorimeter will be in- 
distinguishable from a single-photon shower. To 
estimate the number of QCD background events 
in the Zj sample we use the jets in the Z sample 
and a probability function for jets to fragment 
in such a way that they pass the photon identi- 
fication criteria. This probability function is de- 
termined from events collected by a dijet trigger 
by measuring the rate at which we identify pho- 
tons (as a function of jet Et) and subtracting the 
expected true photon production rate. We also 
correct for the difference in the ratio of jets com- 
ing from quarks versus gluons, since quark jets 
are much more likely to be reconstructed as pho- 
tons. The probability is about 0.29% at 10 GeV 
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FIG. 1. Three-body Mu-y vs. two-body Mu in the 
Z — >■ I-l'^ fJ," and Z — > e^e~ datasets for photons with 
Et>7 GeV. 



and drops rapidly to about 0.07% (with a large 
systematic uncertainty) for jets above 25 GeV. 
The uncertainty on the photon background is 
the dominant uncertainty on the Z"f cross sec- 
tion measurement. 

We find a total of 778 events: 390 e+e"7 
events in 1.07 ± 0.06 fb~^ and 388 fj,~^fi~j events 
in 2.01 ± 0.12 fb-i. We estimate the QCD pho- 
ton background to be 94 ± 26 events (primar- 
ily at low photon Et), and we estimate the 
number of events where a jet is misidentified as 
an electron to be 14 ± 7 events. This gives 
measured cross sections of a{pp — e"'"e~7 -f 
X) = 4.9 ± 0.3 (statistical) ±0.3 (systematic) 
±0.3 (luminosity) pb and cr(pp — >■ ^j.'^ ^~ j + X) = 
4.4±0.3±0.2±0.3 pb. The average cross section 
is a{pp + X) = 4.6 ± 0.2 ± 0.3 ± 0.3 

pb, which agrees with a next-to-leading-order SM 
calculation [17] of 4.5 ± 0.4 pb. These cross sec- 
tions are for events with i?J > 7 GeV, Ai?£-y > 
0.7, and a dilepton mass of at least 40 GeV/c^. 

There is no direct coupling between the Z bo- 
son and the photon in the SM, so Zj production 
at tree level comes from either initial-state ra- 
diation (ISR), where the photon is radiated by 
one of the incoming quarks, or final-state radi- 
ation (FSR), where the photon is radiated off 
one of the charged leptons from the decay of 
the Z boson. In ISR events the dilepton mass 
will approximate the Z mass, whereas in FSR 
events the three-body lepton-lepton-photon mass 
will approximate the Z mass. In Fig. 1 we plot 
the three-body mass on the vertical axis and 
the dilepton mass on the horizontal axis. FSR- 
dominated events form the horizontal line near 
Mwy —91 GcV/c^, and ISR-dominated events 



form the vertical line near M« = 91 GeV/c^. The 
diagonal contains radiative Drell-Yan events. 

Potential anomalous Zj couplings have been 
studied by Baur and Berger [5] and we adopt 
their notation for classifying the nature of the 
coupling: hY where V is either a Z or a 7 and the 
index i runs from 1 to 4. In order to preserve tree- 
level unitarity at large center-of-mass energy -s/i, 
a form factor is used of the form hY — , -';°.2^„ . 
In this paper we use A = 1.2 TeV and n = 3(4) 
for i = 3(4). If an anomalous coupling between 
the Z boson and the photon exists, it will pro- 
duce more events with high-i?J photons than ex- 
pected from the SM. An example of this is in- 
cluded in Fig. 2, which compares the distribution 
of for ISR events (events where Mnj > 100 
GeV/c^) to the SM expectation, showing good 
agreement. 

To set limits on anomalous Zj couplings, we 
use the same dataset used for the cross section 
measurement. Since the statistics are limited 
at high we also include Zj — vDj events 
identified from a photon plus dataset where 
E^ > 90 GeV. The photon plus fJx dataset used 
for this analysis is identical to that used in [ ]. 

The expected signal and backgrounds for 
Z7 — vD^ are described in detail in [] n]. The to- 
tal number of 7 + events with > 90 GeV 
is predicted to be 46.3 ± 3.0, of which 25.6 ± 2.0 
are expected to be Z"f — >■ vvj events; we observe 
40 events. The E'^ distribution is displayed in 
Fig. 3 along with the predicted distributions for 
the signal and backgrounds. 

Using the measured distributions (sepa- 
rated by Z decay channel) we set upper limits 
on the strength of anomalous couplings using a 
binned likelihood method. We determine the ex- 
pected E'^ distributions from Monte Carlo data 
samples with non-zero anomalous couplings as 
well as SM couplings. There are too many sam- 
ples to fully simulate each one, so instead we de- 
termine the net efficiency for a generated event 
to appear in the analysis sample as a function of 
-Ej by fully simulating both SM and represen- 
tative anomalous coupling samples. While the 
difference in efficiency between these samples is 
consistent with statistical fluctuations, we take 
the difference between the measured efficiencies 
as a systematic uncertainty. We then apply this 
efficiency function to generator-level MC samples 
to get the expected E'^ distributions. 

The resulting upper limits on the strength 
of anomalous couplings are shown in Table I. 
The expected limits are determined by generat- 
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ing a large number of pseudo-experiments, Pois- 
son smearing the number of expected signal and 
background events, and varying parameters sub- 
ject to systematic uncertainties (which are dom- 
inated by uncertainties on the background esti- 
mates) [19]. The observed limits are better than 
the average expected limits but fall within the 
expected range. We see no evidence for anoma- 
lous couplings. The limits based upon events 
with Z — >■ e'^e~ and Z — >■ /i^/i^ decays are 
nearly identical to previous limits from the DO 
collaboration [20], which studied the same chan- 
nels using 1 fb~^ of data. Our limits including 
the Z ^ vv decays from 2 fb~^ are not directly 
comparable to the limits published by the DO col- 
laboration [21] using 3.6 fb~^ since the DO limits 
were calculated using A = 1.5 TeV. 

In conclusion, we have measured the Z7 cross 
section at the Tevatron Collider and find that it is 
consistent with SM expectations. We have also 
found that the i?J distribution of photons pro- 
duced in association with Z bosons is consistent 
with SM couplings where there is no direct Z7 
coupling, and we have set limits on anomalous 
gauge couplings. 
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Observed Limit 


Expected Limit 


Limit (/I + e) 


1^3 1 


< 0.050 


0.062 ±0.014 


0.083 


\hl\ 


< 0.0034 


0.0043 ± 0.0009 


0.0047 




< 0.051 


0.064 ±0.014 


0.084 


\K\ 


< 0.0034 


0.0043 ± 0.0009 


0.0047 



TABLE I. Upper limits (95% C.L.) on anomalous Z7 
couplings using notation from reference [">] and A = 
1.2 TeV. The last column gives the limit obtained 
using only the /i^/.i~7 and e"'"e~7 events. 
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